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We present selfdual manifolds for coupled Potts models on the triangular lattice. We exploit two 
different techniques: duality followed by decimation, and mapping to a related loop model. The 
latter technique is found to be superior, and it allows to include three-spin couplings. Starting from 
three coupled models, such couplings are necessary for generating selfdual solutions. A numerical 
study of the case of two coupled models leads to the identification of novel critical points. 

PACS numbers: 05.50.+q, 05.20.-y 

I. INTRODUCTION 

The two-dimensional Potts model is a well-studied model of statistical mechanics 0| and continues to attract the 
interest of many workers. Its definition is simple. Given a lattice with vertices {i} and edges (ij), the Hamiltonian 
reads 

PH^-KY,SiS^:SJ), (1) 

where S is the Kronecker delta. The spins Si = 1,2, ... ,q initially take q discrete values. However, by making a 
random cluster expansion (^j it is easily seen that the partition function can be written 

2 = ^6'=(c)g»(c)^ (2) 
c 

where 6 = e'^ — 1. Here, the sum is over the 2l^*-'^l possible colourings C of the edges (each edge being either coloured 
or uncoloured), e(C) is the number of coloured edges, and n{C) is the number of connected components (clusters) 
formed by the coloured edges. Taking Eq. Q as the definition of the Potts model, it is clear that q can now be 
considered as a real variable, independently of the original spin Hamiltonian. Also, we shall adopt the point of view 
that Eq. 10) makes sense for any real b, although b < —1 would correspond to an unphysical (complex) value of the 
spin coupling K. 

Exact evaluations of Eq. lO, in the sense of the Bethe Ansatz, exist for several lattices and for specific curves in 
space along which the model happens to be integrable This is true, in particular, for the square lattice with 

b - (3) 



and for the triangular lattice with 



b = -2±y/I^, (4) 



b^ + 3b' = q. (5) 



These curves have several features in common. First, they correspond to critical points (with correlation functions 
decaying as power laws) for < g < 4 6], whose nature can be classified using conformal field theory (GET) 0. 
Second, the values of the coupling constants are often so that the partition function is selfdual (see below); this is the 
case for the curves Q and lO above, whereas the two curves in are mutually dual. 

The part of the curves having 6 > corresponds to the ferromagnetic phase transition, whose critical behaviour 
is lattice independent (universal). More interestingly, the antiferromagnetic (— 1 < & < 0) and unphysical regimes 
{b < —1) contain non-generic critical points whose relation to GET has, at least in some cases, not been fully 
elucidated. This is so in particular for b — ~1, where the Potts model reduces to a colouring problem, and Eq. ^ 
becomes the chromatic polynomial. 

Much less is known about several Potts models, coupled through their energy density d{Si,Sj). Results coming 
from integrability seem to be limited to the case of = 2 coupled models which on the square lattice only leads 
to new critical points in the well-studied Ashkin- Teller case (j] (i.e., with q — 2). Apart from that, GFT-related 
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results are essentially confined to perturbative expansions in e ~ <? — 2 around the ferromagnetic critical point |9|, [lOj . 
These results, corroljorated by numerical evidence [TTIIT^. indicate the existence of novel critical points for > 3, 
with possible implications for the random-bond Potts model through the formal analytical continuation (replica limit) 
N ^0. 

In the present publication we investigate the possibility of novel critical behaviour in N coupled Potts models on 
the triangular lattice. To identify candidate critical points we first search for selfdual theories. In comparison with a 
similar study on the square lattice 0,0] several distinctive features emerge due to the non-selfdualness of the lattice. 
This leads us to use two different techniques. In the first, a standard duality transformation is followed by decimation 
(star-triangle transformation). This turns out to be quite cumbersome, already for N — 2. We therefore turn to a 
second technique, which utilises a mapping to a system of coupled loop models. This leads to simpler relations, and 
as a bonus allows to include three-spin couplings around one half of the lattice faces. Starting from = 3 coupled 
models, such additional couplings are actually necessary for generating non-trivial selfdual solutions. 

For TV = 2 we numerically investigate the non-trivial selfdual manifold. Following the motivation given above, the 
main interest here is to establish whether a given selfdual point corresponds to a renormalisation group fixed point 
(and possibly even to a critical fixed point). We shall see that these expectations are indeed born out: the numerics 
is compatible with critical points whenever < g < 4. Measuring the central charge, we identify the corresponding 
universality classes. These can in some cases be understood from those of a single model, but we also identify points 
possessing novel critical behaviour. 

The paper is laid out as follows. In Section^ we present the technique of duality followed by decimation for two 
coupled models with pure two-spin interactions. In particular, we find a non-trivial selfdual solution. The mapping 
to a loop model, given in Section [llll allows to rederive this solution in a much simpler way, and to generalise to the 
case where three-spin interactions are included. In Section Hvl we use this technique to treat the case of three coupled 
models with both two and three-spin interactions. A numerical study of the non-trivial selfdual solution found in 
Section is the object of Section Finally, Section IVll is devoted to our conclusions. 



II. MODELS WITH TWO-SPIN INTERACTIONS 



To illustrate the first technique (duality and decimation) , we consider the case of iV = 2 coupled models with two- 
spin interactions. In order to simplify the notation, wc introduce the symbol (5^ = S (S'f , Sj^, where the superscript 
refers to the spins of the /i'th model (p = 1,2,...,N). We are interested in the coupled model defined by the 



Hamiltonian 



0H2 = -J2 {KiSl, + K^Sl + K,25},5l } . (6) 

(ij) 

The spins S'f take different values. 

A. Duality followed by decimation 

As shown in Ref. [l^ Il3| , Eq. ((BJ admits a (generalised) random cluster expansion resulting in 

Z=Y, ^e(C,nCj)^e(CrnC.)^e(C,nC.)^«(C.)^«(C.)^ 

Cl ,C2 

where are independent colourings of the /i'th model, and we have defined the complementary colouring = 
{ij) — Cfj_. The new parameters b are related to the coupling constants K through 

= e^- - 1, bi2 - c^'i+^"^+^i^ - e^i - e^^ + 1 (8) 

As explained in the Introduction, we shall take the point of view that the model is defined by Eq. ((7|) for any real 
values of b and (7^. 

Up to an irrelevant constant, the partition function of the dual model is again given by Q, but now with respect 
to the dual (hexagonal) lattice, and with dual values b of the parameters ;12.,13i|: 

-6291 r biq2 ~ 5192 

^1 = -^, &2 = ^, 612 = ^. (9) 




FIG. 1: The star-triangle transformation. 



A rather obvious procedure would be to follow @ by a standard decimation prescription (star-triangle transformation) 
in order to get back to parameters b' defined with respect to the triangular lattice, and then search for selfdual solutions, 
b = b' . A key assumption, of course, is that such solutions exist within the original parameter space, i.e., with only 
nearest- neighbour couplings among the spins (141 . 

The precise setup is shown in Fig. ^ We form the partial trace over all spins Sq situated at even (F-shaped) 
vertices of the hexagonal lattice, while keeping the exterior spins S'^, S2, fixed. Defining b as in Eq. (jSJ, we must 
have 

^ exp 1^ (k, ^ Si;, + K;2S'o.S'o)j | = Aexp J ^ ( i^^ ^ 5^^ + K[,Sl^6l) I , (10) 

where the proportionality factor A does not have any bearing on the duality relations for the coupling constants. 

We obtain ten distinct relations by considering all symmetry-unrelated choices for the fixed spins S'f with fj, — 1,2 
and i = 1, 2, 3. Following 14] we suppose q^, > 3 integer initially, and then invoke analytic continuation to claim the 
validity of the result for arbitrary q^. 

• For 5*1 7^ 5*2 7^ 53 on the two lattices: 





-3)(g2-3) + 3(l + 6i+&2 + foi2)- 


1-3(92 


-3)(l + 6i) + 


3(gi- 


-3)(l + 52) 


+ 6(1 + 


bi){l+b2)=A 


(11) 


• For 5*1 


= S2 = S3 on the two lattices: 
















(91 


- I)(g2 - 1) + (1 + 61 + ;r2 -1- b^2f - 


f (92- 


-l){l + bif + 


(91- 


m + b2f 


= A(l + 


-b[+b'2 + b[,f 


(12) 


• For Si 


= S2 ^ S3 on the two lattices: 
















(91 


- 2){q2 - 2) + (1 + bi + b2 + b^2f 


+ (1^ 


- 61 + 62 + bi2] 


+ (92 




' + (91 


-2)(l + 6~2)2 + 






{q2-2){\ + bi) + {qi-2){\ + 


62)(1 


f 6i)'(l + 6~2) 


+ (1^ 


-62)^(1 +61 


)=A{1 


+ b[ + b',+b[,) 


(13) 



• For Si = S^ ^ Si and 5? ^ S^ = S^: 

(91 - 2)(<72 - 2) + {q2 - 2)(1 + hf + (92 - 2)(1 + h) + (91 - 2)(1 + h)^ + [qi - 2)(1 + h) + 
(1 + bi + h + &r2)(l + 6'i)(l + 62) + (1 + bi + h + &l2)(l + bi) + (1 + bi){l + h) + 

(1 + h + h + h2){l + h) = A{1 + b[){l + b'^) (14) 

• For Si ^ S^ = Si and Sf ^ S^ ^ S^. 

(91 - 1)(92 - 3) + (92 - 3)(1 + + 3(91 - 1)(1 + h) + 3(1 + bi + h + &r2)(l + bi)^ = A{1 + b[f (15) 

• For Si = 51 = 51 and Sf = Si ^ Sf. 

(92 - 2)(9i - 1) + (91 - 1)(1 + b2f + (91 - 1)(1 + b2) + (92 - 2)(1 + bif + 
(1 + 6~i + &~2 + b^2f{l + 6i) + (1 + bi + 6~2 + &T2)(1 + bif = A{1 + b\ + b'^ + 6'i2)(l + ^i)" (16) 

• For Si = Si ^ Si and 5? ^ Sj ^ Si- 

(91 - 2)(92 - 3) + (92 - 3)(1 + 6~i)2 + (92 - 3)(1 + 6~i) + (91 - 2)3(1 + 61) + 
2(1 + 6~i + &~2 + &r2)(l + bi) + (1 + &i)2(l + 62) + (1 + 6i + b2 + ^2) + 2(1 + bi){l + 6~2) = ^(1 + fo'i) (17) 

• Eqs. (|15|l . Hlfci|) and H17|) each represent a pair of relations of which we have only written one representative; the 
other one is obtained by permuting the two models, i.e., by letting 91 ^ 92 and bi ^62- 
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FIG. 2: Variation of the Boltzmann weights l|2()^ with the parameter g defined in Eq. H19^ . 

B. General structure of trivial solutions 

The list of selfdual solutions of N coupled models, fi — 1,2, ... ,7V, will in general contain a certain number of 
trivial solutions. By a trivial solution we mean one which is a consequence of the selfduality of a single model. Let 
us discuss in detail two classes |0, 0| of trivial solutions: 

1. The models are actually decoupled. This happens, e.g., in the above example with N = 2 when K12 = 0, that 
is 5i2 = 6162- We will then have 6^ + 36^ — for each /i, cf. Eq. ||SJ). The number of real solutions of the ^'th 



equation is = 3 when < < 4, 



be n = njy^]^ trivial solutions for the system of N decoupled models. 

2. The models couple strongly so as to form a single q — n^=i 9^ state model. This happens when only the 
coupling constant involving all the models is nonzero. E.g., in the above example with iV = 2, one would have 
Ki — K2 = 0, that is 61 = 62 = 0. The number of such solutions equals the number of real solutions of Eq. lO, 
with b replaced by 612. 

The goal of our study is to show that there exists selfdual solutions of coupled Potts models on the triangular lattice 
which are not trivial in the above sense. 



2 when (7^ = or = 4, and = 1 otherwise; there will therefore 



Non-trivial solutions 



Let us return to the Hamiltonian jnj. We have numerically solved the ten relations (|ll(l - (|17|l for several different 
values of qi and 92- The conclusion is that for qi ^ 92 only trivial solutions exist. 

For qi — q2 the situation is different. There are now only seven distinct relations (|11|I - H17() . since the three relations 
which formerly occurred in pairs will now collapse into single relations. The parameters b are thus less constrained, 
and accordingly we find non-trivial solutions. Numerically we find that these solutions have 61 = 62 (but note that 
there are still trivial solutions with K12 = which break this symmetry). 

Setting now q = qi = q2 and 61 — 62 we can obtain the non-trivial solutions analytically, e.g., by solving Eqs. lll|l . 
(|13|1 and IjlTf) for A, bi and 612, and verifying that the found solution satisfies all the other relations. The result is: 



bl + 6bl + 3biq + q{q - 2) = 0, 



(18) 



For each q e (0,4) Eq. H18|) admits three distinct solution for 61. To make clear in the following exactly to which 
solution we are referring, it is convenient to recast (|18|l in parametric form, by setting q = Acos^{-ng). When the 
parameter g runs through the interval [0, |], the number of states q runs through the interval [0,4] three times. We 
have then 



bl ~ x{l — x), 612 — {x — 1)^ 



o /27r 
2 cos ( —g 



(19) 
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This parametrisation also has the advantage over (P|l that it is non-singular as 5 ^ 1 (i.e., x — > — 1 and q ^ 4) and 
yields the correct limiting values of bi and 612 . 

In terms of x, the Boltzmann weights for two neighbouring spins being identical in none, one, or both of the two 
models read: 

1, e-^i = 1 + a; - a;^ e^^^+^^i = 2 - x^. (20) 

Their variation with g is shown in Fig. |2] Note that the K12 coupling is physical (e^^^ > 0) for | < <? < |, and that 
Ki is physical for ^ < g < ^. 

It is also interesting to remark that in the x-parametrisation, Eq. (|3J| for a single model reads b ^ x — I; the trivial 
solution of type 1 is then 61 = a; — 1, 612 = (61)^ with the same value of 612 as in Eq. (fT^ . 



D. Special points on the curve (1191) 

Let us remark on a few special values of the parameter g for which the physics of the two coupled models can be 
related to that of a single model. 

1. For g = 1 (i.e., q = 4 and x = —1) one has K12 = 0, whence the models are decoupled. 

2. For .9=1 (i.e., q — 2 and a: = 0) one has Ki — and K12 — log 2. Thus, the two models couple strongly to 
form a single = 4 state model at the ferromagnetic fixed point. 

3. For g = ^ (i-c, q = and x = 1) one has Ki = K12 = 0. This is an infinite-temperature limit, whose properties 
depend on the ratio K1/K12 a.s g ^ ^. In fact, for x ^ 1 we find Ki ~ (1 — a:) and K12 ~ 2(1 — x)^, whence the 
coupling between the two models is negligible. Note also that q ^ 3(1 — a;)^, whence the point (g, 61) = (0, 0) is 
approached with infinite slope, as is the case for a single Potts model along the selfdual curve ISJ. 

Note also that these values of g correspond to degeneracies of the Boltzmann weights H2()(l , cf. Fig. [21 

As to the critical behaviour of a single Potts model, the situation is well understood in the case of the square lattice 
0. There are three critical phases, referred to as ferromagnetic, Berker-Kadanoff and antiferromagnetic in Ref. 0. 
By universality, one would expect the same three critical pha ses to describe the distinct branches of the selfdual curve 



pnas 

as is confirmed by numerical transfer matrix results |l6| . In particular, for the central charge along jSJ one has 

63^ 1 

for < g < — (ferromagnetic phase) 



1-9 --2 

1 

c = 1 , for - < g < 1 (Berker-Kadanoff phase) 

1- g 2 ■ 

3 

c — 2~ 6{g — 1), for 1 < 5 < - (antiferromagnetic phase) (21) 

when parametrising g = 4cos^(7r5) as in H19|l . Note that points (g, 61) = (0,0) and (g, foi) = (4,-2) are special, and 
the critical behaviour when approaching the curve (O at these points depends on the exact prescription for taking 
the limit. 

We conclude that the above-mentioned special points on the selfdual curve p9|l for two coupled models lead to the 
central charges c = —4 for g — ^ and c ~ 1 for .g = |. For 5 — > 1, the result is for the moment uncertain, due to the 
ambiguity in taking the limit just referred to. 



III. TWO AND THREE-SPIN INTERACTIONS 



The results of the preceding section can be generalised to a model defined by the Hamiltonian H = H2 + H3. The 
term H2 is as in Eq. ©, whereas introduces interactions between the three spins around the up-pointing faces 
{ijk) of the triangular lattice, as shown in Fig. 01 

- ^ + L24, + L,2Sl,kSf,k } • (22) 

(ijk) 

Introducing such interactions around every face seems difficult, even in the case of a single model 0. 
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FIG. 3: The lattice of Potts spins is shown in solid linestyle. There are two-spin interactions along every edge, and three-spin 
interactions among the spins surrounding the up-pointing faces (shaded). The loop model is defined on a shifted triangular 
lattice, shown in dashed linestyle. 

One could of course consider also interactions between two spins in one model and three spins in the other. It can 
be verified that the method given below can be adapted to this case. However, we have chosen not to consider any 
further such mixed interactions. 

In the following it is convenient to introduce the parameters 

y^ = e^^-l, yi2 ^ e^''+^'+^' - e^' - e^^ + 1 (23) 

in analogy with Eq. ||SJ|. Guided by our results without three-spin interactions, we shall assume in the following that 
non-trivial selfdual solutions only exist when coupling identical models. Thus we restrict the study to q = qi =92, 
61 — 62 and Li ^ L^. 

A. Mapping to a fully-packed loop model 

Wu and Lin IT5| have shown how to produce duality relations for a single Potts model with two- and three-spin 
interactions, by mapping it to a related loop model. After briefly reviewing their method, we shall adapt it to the 
case of coupled models. 

In Fig. Owe show the triangular lattice of Potts spins, and the shifted triangular lattice on which the loop model is 
defined. To obtain the correspondence, one first rewrites the Boltzmann weight around an up-pointing triangle (ijk) 
as 

5 

a=l 

where Si = 5ij ~ 5{Si, Sj) etc. Note that =^ 1. To each of the five terms in this sum is associated a link diagram 
on (ijk), as shown in the first line of Fig. 01 indicating which spins participate in the delta symbol. The partition 
function is then 

^potts = E n -'.^ - E n /^^'^ (25) 

{i]k) C a=l 

where the sum is over all link diagrams C for the whole lattice, n(£) is the number of connected components in £, 
and na{C) is the number of up-pointing triangles whose link diagram is of type a = 1, 2, . . . , 5. 

The link diagrams are now mapped to fully-packed loop configurations on a shifted triangular lattice (cf. Fig. ^ 
via the correspondence given in the second line of Fig. ^ The partition function of the loop model is defined as 

^ioop=E-^^''^ric^'^ (26) 

a a=l 

where p(£') is the number of closed polygons (loops) in the loop configuration £' (which is in one-to-one correspondence 
with C using Fig.^, and z is the fugacity of a polygon. Using now the Euler relation one finds that ^potts = ^loop if 
z ^ y/q, Ca = fa for a = 1, 2, 3, C4 = 9^/^/4, and C5 = q'^^^fz El- 

Finally, the duality relation follows from the invariance of the loop model under 7r/3 rotations; notice that this 
cyclically interchanges the link diagrams of types 1, 2, 3 and permutes the diagrams of type 4, 5. 




FIG. 4: Correspondence between link diagrams L for the Potts model (first line) and vertices £' of the fully-packed loop model 
(second line). 

B. N = 2 and mapping to coupled loop models 

The above mapping can be adapted to the case of N coupled models. We here consider N — 2. The Boltzmann 
weight Wijk around an up-pointing triangle (ijk) can be written in a form generalising H24() : 



Wi-jk 



2 . 

b 1 



(27) 



,,b=l 



as usual (5'' refers to the /i'th model. As the two models are identical, we have the symmetry fab — fba- When 
Li = L2 = L12 = 0, the couplings are denoted /"j,; they are related to the b through 



J i4 
fO 
Jib 

J 45 

nO 

J55 



= "12, 

= bl, 



bi2{bi + 2h) + bl 
1, 

bl + 3bl, 

6?2 + 36?2(l + 2&i) 



h2bl 



(28) 



where i j are any numbers in {1, 2, 3}. In the general case (L 7^ 0) we then have 

/l4 = /i4, 

/l5 - /l"5+yi(/?i + 2/i"2 + /l"4 + /?5), 



/44 — /° 



44, 



/55 = /5°5 + 2yi(3/°5 + fl + /°5) + yi2(3/°i + e/f^ + 6A"4 + Qfl + /^^ + 2fl + f,,) : 

note that we here only give the fab needed in the duality relation H33|) below. 
The Potts model partition function reads 



(29) 



^Potts 



n(Ci)+n(C2) 



n/ 

a, 6=1 



n„b(£i,£2) 
ah 



(30) 



Here, nab{Ci, C2) is the number of up-triangles where model 1 is in the link state a and model 2 is in the link state 
b. On the other hand, we can define a coupled loop model through 



n 

a, 6=1 



n„b(£i,£2) 



(31) 



The equivalence between the two goes through as before. Using the Euler relation, one finds .^potts = ^loop provided 
that z = y/q and that 



„ , _ „(l5a,4+<5b,4-<5o,5-<56,5)/2 f , 

t^ao — H J ab- 



(32) 
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It should now be obvious how the mapping generalises to an arbitrary number of coupled models. One has Cat- - = 
^(Af4-A'5)/2y^^ ^ where -/V4 (resp. N^) counts the number of indices equal to 4 (resp. 5). 

The selfduality criterion is again obtained by requiring invariance of .^loop under 7r/3 rotations. This means that the 
Cab are invariant under independent permutations of the indices {1, 2, 3} and of {4, 5}. We also recall the invariance 
under a permutation of the two models, Cat — Cha- Actually, since the two Potts model were taken to be identical and 
isotropic from the outset, the only non-trivial selfduality criteria are: 

C44 = C55, ci4 = C15. (33) 



C. Selfdual solutions 



We now wish to express the condition of selfduality in terms of the parameters h and y, cf. Eqs. ^ and H23|l. 

When three-spin interactions are absent [L — Q), the relations H33() immediately yield the solutions given in Sec. Ill Cl 
This is a remarkable simplification when compared to solving the system of relations Hll(l - (|17|l : indeed many of these 
relations turn out to be dependent. 

When L 7^ 0, Eq. Ip!^ still gives the complete solution to the selfduality problem, but it does not generically lead 
to simple expressions in terms of the parameters b and y. We therefore concentrate on a few remarkable solutions. 
As before, there are two types of trivial solutions: 

1. Trivial decoupled solution {K12 — £12 = 0, or &12 = b\ and yi2 — yf). One finds the selfduality criterion of a 
single model 

bl + 3bj-q + yiil + bif = 0. (34) 

2. Trivial strongly coupled solution {Ki = Li = 0, or bi = yi = 0). We find 

5?2 + 36?2-'?'+yi2(l + 6i2)' = 0. (35) 
This is just the selfduality criterion of a single state model. 
Some noteworthy non-trivial solutions can be found by giving particular values to 5i, yi or j/12: 

1. For yi = (i.e., Li — 0), there is only one non-trivial solution: 

{bl + 6bj + 3qbi + q{q-2)){bl + ^bl-q) =yi2ibi + 5bi + q + 2)^, 612 = (36) 

Note that when j/12 — 0, the factorisation of the left-hand side allows us to retrieve either (|18ll or the trivial 
solution 

2. For yi2 — yf (i.e., L12 = 0), there is a solution of the form: 

h 1 h 9^ , Q Q 9(4 -g) ,„„^ 



3. For 61 = —1 (i.e., Ki — > —00), there is a solution of the form 

2 (^^■^ ■ q^ -Aq + 4 



1 / q'-5q + 5 \ 

yi = -7: yi2 + — ] — — ■ (38) 



IV. THREE COUPLED MODELS 



The technique of mapping to coupled loop models has permitted us to study the case of = 3 coupled Potts 
models, defined by the Hamiltonian 



E W E + E + Li235l,u5lu6%A . (39) 

{ijk) I ^=1 ii>u=l ) 
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Since in the case of two coupled models, nontrivial selfdual solutions were only found when coupling identical models 
in an isotropic way, we shall henceforth restrict the coupling constants as follows: 

Ki=K2 = K3, Ki2 ^ Ki3 ^ K23, Li=L2 = L3, il2 = ii3 = -^23, (40) 

and we will use the parameters 0, 0| 

61 = e^^ - 1, 612 = e^^^+2^^ - 2e^i + 1, 6123 = e^-^+^^^+^^i - 36^^^+^^^ + Sc^^ - 1 (41) 
2/1=6-^1-1, yi2 = - 2e^i + 1, 2/123 = e^i"^+''^^i"+'''^i - Se^^'+^^i + 3e^i - 1. (42) 

The mapping to coupled loop models follows the obvious generalisation of Eqs. (|30|I - H31|I . and the equivalence 
criterion is stated in the remark after Eq. H32|l . 

To relate the coupling constants fabc to the b and y, we first consider the case of vanishing three-spin interactions 
(i.e., y = yi2 = j/123 = 0). Letting i ^ j ^ k designate distinct numbers in {1, 2, 3} we have: 



J III 


— bi23, 






= ^1^12, 




rO 

Jijk 


- bl 




fO 


= bi2, 




M 
•> ii5 


= bi23b\ 


+ 2612361 +5l25?, 


fO 

Jij4 


= bl 




rO 

Jij5 


= bUl 


+ 61) + 26126?, 


fO 


= bl, 




rQ 

J lib 


= buibj 


+ 26i) + 6?. 


M 

Jihb 


= 25?2 4 


- 56^361 + 26123612 + 46i23' 


rO 

J444 


1, 




rO 

J445 


= 5? + 35?, 


J455 




3622 + 661262 + 662261, 


fO 

J 555 


= 65?2 4 


- 18612(61 + 612)6123 + 3(1 



(43) 



61261 + 261236? + 61236^2, 



36i+ 3612)6^23 + ^>l23- 

For the general case one then has (note that we only give those fabc which will be used in the duality relation H45f) 
below) : 

/ll4 = /1147 

/ll5 = /?i5 + yi(/?il+2/?i2 + /?i4 + /l"l5), 
/l24 = /1247 

/l25 = /l25 + yi(2/ll2 + /l23 + /l24 + /125)' 
/l44 — /i44, 

/l55 = /?55 + Mfll5 + 2/?25 + fl45 + A^Ss) + ^12 (/ni + /?44 + fl55 + 6/{'l2 + 2/?23 + (44) 

2/{'l4 + 2/i°i5 + 4/i"24 + 4/025 + 2/1O45), 
/444 = /444 7 

/445 — /445 + y(3,/l44 + /444 + /445)' 

/455 = /4"55 + 2y(3/°45 + /4°45 + fD + 2/l2(3/?i4 + ,/4°44 + /4°55 + 6/i°24 + 6/i"44 + 6/i°45 + 2/4°45): 

/555 = /5"55 + 3y(3/i°55 + fl, + ^"gg) + 'Syui^^i^ + f445 + /555 + 6/?25 + 6/i°45 + ^.fl55 + V455) + yM^flll + 
/444 + /555 + 18/?12 + 9/fi4 + 9/i°i5 + 6/1023 + 18/iV + 18/i°25 + 9/?44 + 9/i°55 + Wl45 + 3/4°45 + 3/4°55)- 

So from (the generalisation of) Eq. H32|l we know how to express the Catc in terms of the 6 and y. Given that we 
have isotropic Potts models with the same coupling constants, the non-trivial selfduality relations read simply: 

C444 = C555, C455 = C445, C155 — C144, C115 = C114, C125 = C124. (45) 
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A. Selfdual solutions 

There are two types of trivial solutions, as in the case of two coupled models. 

An important difference with the case of two coupled models is that non-trivial selfdual solutions with y = j/12 = 
2/123 = only exist for exceptional values of q (i.e., g = and (7 = 4; see below). So for generic values of q, the 
three-spin interactions are necessary to generate non-trivial solutions of the selfduality problem. 

As before, the most general solutions are not algebraically simple in terms of the variables h and y. We therefore 
report only a few special cases: 

1. There are three non-trivial solutions with L12 = ^123 = (i.e., yi2 = yl and 1/123 = Ui)- They read: 

, 9 , , - 9g' + i8g - 12 g(4 - g) 

h = --. 612 = ^, 6123 = ^ , y^- q^-Aq + A '^ 

61 = -1, 612 = |, ^1.3 = ^, yr^^; (47) 
.2 , , „_n . _ '7(l + 2&i) 



(4g - 6)6f + 2g6i + g = 0, 612 = -- 



2 

g^((8g^-16g-6)bi+4g^-12g + 3) q 
4(36i+,)(2,-3) ' 2—q 

2. There are two non-trivial solutions with Li = Q (i.e., yi =0): 

fci = -1, ^ |, &123 = '^^^ ^\ V12 = -r—-r^^ yi23 = 2 ,^!o^^''n ; (49) 

2 2 — 4g -|- 4 gr-^ — bq'^ + 12q — 8 

, , , g(l-g) , q{q^ - 3g + 1) 

01 = 1-9, 012 = ^^ , 0123 = 7:^ TT ^, 

2-(7 (2-(7)((7-3) 



yi2 



96 - I8g5 + 126(?4 - 432^3 + 756g2 - 648(j 216' 



_ g(4gS - 72g^ + 531g^ - 2068g^ -I- 4584g^ - 5856g'^ -I- 4220g^ - 1584g + 240) 

^^^^ ~ g9 - 30g8 -I- 390g7 - 2872g6 + 13140^5 - 38520g4 + 71928g3 - 82080g2 + 51840g - 13824' ^ ' 

3. There are two non-trivial solutions with yi2 — 0. We give here only the first one because the other is complicated: 

, , , 9 , g(l - g) g(4 - q) q\q-&) 

= ^^^-2' ^123 = ^—, 2/i-2(g2_4<z + 4)' ^^^^-2(g3_ 69^ + 129 -8)- ^^^^ 

Note that for q = 2 the non-trivial solutions given are singular. In fact, for q ~ 2, they can be written as: bi — — 1, 
612 = 1, 6123 — —1, and the values of y, yi2, ?;i23 are arbitrary. Indeed the values of the b correspond to three 
decoupled antiferromagnetic Potts models at zero temperature, and so the values of the parameters yi, 2/12, 2/123 do 
not matter. 



V. NUMERICAL STUDY OF TWO COUPLED MODELS 



It was mentioned in the Introduction that the Potts model is usually critical on the selfdual manifolds, for suitable 
values of q (i.e., < 9 < 4). We expect this also to be true for coupled Potts models, and so it is interesting to 
determine the corresponding universality classes. In the lack of an exact (Bethe Ansatz) solution, this question can 
be addressed by evaluating the effective central charge along the selfdual manifolds, e.g., using numerical transfer 
matrix techniques. 

In this Section we focus on the selfdual curve p9|) for two coupled models with pure two-spin interactions. Note 
that in Section [ll Dl we have already remarked on a few special values of the parameter g for which the physics of the 
two coupled models can be related to that of a single model. 
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FIG. 5: Semi-infinite strip, here of size L = 2 triangles in the finite (vertical) direction. Periodic boundary conditions identify 
the top and the bottom of the figure. Potts spins are defined at the loci of the small circles; they interact along the solid lines 
which form a triangular lattice. The labels within each circle identify the usual three sublattices of the triangular lattice. The 
loop model is defined on the medial Kagome lattice, shown in broken linestyle. The transfer matrix propagates the system 
along the horizontal direction, from left to right. Thin dotted lines indicate successive time slices (see text). 



A. Transfer matrix 



The triangular-lattice Potts model can be transformed into a loop model on the medial (surrounding) graph — which 
is the Kagome lattice — in a standard way JlJ. (This loop model should not be confused with the one described in 
Section Unl l We have computed the free energy of the two coupled models © on semi-infinite strips by constructing 
the transfer matrix of two coupled Kagome-lattice loop models. 

The geometry is depicted in Fig. [S] For a periodic strip of circumference L triangles, each time slice cuts 2L 
dangling ends of the Kagome-lattice loop model. In order to have the leading eigenvalue Aq of the transfer matrix 
Tl correspond to the ground state of the continuum model, the definition of Tl must respect the usual sublattice 
structure of the triangular lattice. This means that L must be even, and that successive time slices are as shown in 
Fig.m 

The numerical diagonalisation is most efficiently performed by decomposing in a product of sparse matrices, 
each adding one vertex of the Kagome lattice (or, equivalently, one edge of the triangular lattice). With the setup of 
Fig. |31 all these sparse matrices are identical, except for the position of the two dangling ends on which they act. We 
have been able to diagonalise Tl for sizes up to L = 10 (the corresponding matrix has dimension 141 061 206). 



B. Central charge 

The free energy per unit area is f{L) = — ^J^j^ log^Oj the length scale being the height of one triangle. We have 
extracted values of the effective central charge c from three-point fits of the form ^3 

/W = /M-^ + ;^> (52) 

where the non-universal term in A is supposed to adequately represent the higher-order corrections. 

Fig. El shows the numerical values of c{g) along the curve H19(l . For each value of g, three estimates for c{g) are 
shown, obtained by fitting {/(L - 4), f{L - 2), /(L)} to for L = 6, L = 8 and L = 10 respectively. 

Naively, one would expect the K12 coupling to be marginal at the q — 2 ferromagnetic point (i.e., at g ~ |) 
and the surrounding regime to be accessible to perturbative calculations. However, it should be remembered that 
1) the point g — j is not that of two decoupled Ising models but that of a single 4-state model, and that 2) the 
renormalisation group equations for N coupled models are singular when N = 2 Nevertheless, the numerics 

seems quite conclusive that the q < 2 regime with j < g < j has a central charge which is just twice that of 1)21(1 
[upon changing the parametrisation, g ^ 1 — <?] : 

c(.9) = 2 (1 - ^ii^) , forl<,<J, (53) 



meaning presumably that the continuum limit is really that of two decoupled models. This is also consistent with 




FIG. 6: Three-point fits for tfie effective central cliarge c as a function of g. Tlie soiid curve sliows tfie exact resuit 1)53^ . valid 
for I < 5 < |, as discussed in tfie text. 



the result c{g = \) ~ —4. The agreenient of the numerics with (|53ll is excehent also for those data points (with 
0.25 < g < 0.35) which are not visible in Fig. 

The region < g < gi, with gi « 0.15, is interesting as it corresponds to c > 2. This hints at the coupled 
models requiring some kind of higher symmetry than its two constituent bosonic theories. Note in particular that 
for g = 0, our three estimates for c read C2,4.6 = 2.758, c^^e.s — 2.914 and cg.s.io = 2.966, which we extrapolate to 
c{g = 0) = 3.00 ± 0.01. We conjecture that the exact result is c(a = 0) = 3. Since q ~ i, this theory can also be 
represented as two coupled vertex models on the Kagome lattice Q • 

For g = 1, the Boltzmann weights H20I) are all ±1. It turns out that in this particular case it is more convenient to 
work with a modified transfer matrix that adds not one but L/2 time slices, cf. Fig. |S1 This matrix has its largest 
eigenvalue equal to unity regardless of L, and we conclude that f{L) = for any L. In particular, this means 
c{g = 1) = 0. The numerics is however indicative of a non-trivial regime for | < g < 1, and it seems that we may 
have c{g) — > 4 as g — > 1~, consistent with two decoupled models each of which is obtained by taking the limit 5 — *■ 1"*" 
in lEH). 

In the region I < g < g2, with 52 ~ 1.10, our numerical diagonalisation scheme experiences difficulties, maybe due 
to the leading eigenvalue having a non-zero imaginary part. 

Finally, in the regime 52 < <? < | the central charge takes large, negative values (in particular, some of the values 
are not visible in Fig.EJ. At first sight one might believe that the continuum limit is that of two decoupled models in 
the Berker-Kadanoff phase, i.e., with c{g) given by twice that of (|21|l [upon changing the parametrisation g 2 — g]: 
c{g) = 2(1 — 6(2 — gY /{g — 1)). This possibility is however clearly ruled out by the numerics, and c{g) appears to be 
given by a non-trivial expression. 

As g — > |, the two leading eigenvalues of the transfer matrix become degenerate. In the sense of analytically 
continuing the curve (|18|l to negative values of q, this presumably marks a transition to non-critical behaviour for 
q < 0, i.e., with the phase transitions being first-order in b upon crossing the curve p8(l . 
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VI. CONCLUSION 



Using a mapping of coupled Potts models on the triangular lattice to coupled loop models, we have obtained non- 
trivial selfdual manifolds for two and three coupled Potts models with two and three-spin interactions. A numerical 
study of the case of two coupled models shows that these manifolds are good candidates for novel critical points, in 
particular in the antiferromagnetic and unphysical regimes. 

The technique can be applied to any number of coupled models, but expressing the solutions explicitly in terms of 
the original coupling constants becomes increasingly complicated as the number of models grows. This is in contrast 
to the quite simple results for coupled Potts models on the square lattice 0, 0| . 

It would be interesting to study the simplest non-trivial case (|19|l using the methods of integrable systems. 
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